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Abstract 

In this paper we consider two semimartingales driven by diffusions and 
jumps. We allow both for finite activity and for infinite activity jump compo- 
nents. Given discrete observations we disentangle the integrated covariation 
(the covariation between the two diffusion parts, indicated by IC) from the 
co-jumps. This has important applications to multiple assets price modeling 
for forecasting, option pricing, risk and credit risk management. 
An approach commonly used to estimate IC is to take the sum of the cross 
products of the two processes increments; however this estimator can be highly 
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biased in the presence of jumps, since it approaches the quadratic covariation, 
which contains also the co-jumps. Our estimator of IC is based on a threshold 
(or truncation) technique allowing to isolate all the jumps in the finite activity 
case and the jumps over the threshold in the infinite activity case. We prove 
that the estimator is consistent in both cases as the number of observations 
increases to infinity. Further, in presence of only finite activity jumps 1) IC is 
also asymptotically Gaussian; 2) a joint CLT for IC and threshold estimators 
of the integrated variances of the single processes allows to reach consistent 
and asymptotically Gaussian estimators even of the /3s and of the correlation 
coefficient among the diffusion parts of the two processes, allowing a better 
measurement of their dependence; 3) thresholding gives an estimate of IC 
which is robust to the asynchronicity of the observations. 
We conduct a simulation study to check that the application of our technique 
is in fact informative for values of the step between the observations large 
enough to avoid the typical problems arising in presence of microstructure 
noises in the data, and to asses the choice of the threshold parameters. 

Keywords: co-jumps, integrated covariation, integrated variance, finite 
activity jumps, infinite activity jumps, threshold estimator. 

1 Introduction 

We consider two state variables evolving as follows 

dX { t q) = a {q) dt + ai q) dW t (q) + dJ^\ q=l,2 
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for t G [0, T] , T < oo fixed, where a and a are cadlag stochastic processes; Wt = 
PtWf 1 ^ + \/l — p% Wf 3 ^; and are independent standard Brownian motions; 
and and are possibly correlated pure jump semimartingales. Given discrete 
observations X^ ,Xof, with observation times spanned on [0, T], we are interested 
in the separate identification of the integrated covariation ICt '■= Jq ptcr^a^dt, 
between the two diffusion parts, and of the co-jumps AJ t A J t , the simultaneous 
jumps of X^ and X^ 2 \ where, for each q — 1,2, AJ^ denotes the size J t — J t _ 
of the jump occurred at time t. 

The recent empirical interest on co-jumps in financial econometrics is motivated 
by the problem of a correct assets price model selection. This has important con- 
sequences in forecasting, in option pricing, in portfolio risk management, and even 
in the credit risk management, since a default of a firm is interpretable as a jump 
in the firm value and contemporaneous defaults give a co-jump, implying default 
dependence (contagion, [1 0J ) . 

ijp (\\ ('2'] 

A commonly used approach to estimate j Q p t cy\ o~\ dt is to take synchronous 
and evenly-spaced observations X^\ X^\ ...X^\ XP\ X^\ ...X^\ with t n = T, 
and to consider the sum of cross products YJ l j= i^jX {l) ^jX^\ where AjX® : = 
x[^ — Xj;^; however this estimate can be highly biased when the processes X^ 
contain jumps; in fact, as n — > oo, such a sum approaches the global quadratic 
covariation 

[X«X (2) ]t= I* Pt *l 1) a? ) dt+ £ AJ«Ajf , 

•*° 0<t<T 

which contains also the co-jumps. To our aim it is crucial to single out the time 
intervals where the jumps occurred. 

A jump process J is said to have finite activity (FA) when a.s. only a finite 
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number of jumps can occur in each finite time interval. On the contrary J is said 
to have infinite activity (IA). In the special case where J is Levy and has IA then 
a.s. infinitely many jumps occur in each finite time interval. 

Our estimator of ICt is based on a threshold criterion allowing to identify all 
the time intervals }tj-i,tj] where the path of a univariate semimartingale jumped, 
if the jump component J has FA, and the intervals where jumps over the threshold 
occurred, if the discretely observed realization of J has infinite activity. Extending 
the application of the criterion to a bivariate framework allows to derive an asymp- 
totically unbiased estimator of ICt as well as of the co-jumps occurred up to time 
T. More precisely we construct the following estimator 

n 

IC T , n '■= ^jA (1) l| (A ^ X (i) )2 < r(ft) }A : ,-A (2) l|( Aj . X (2))2< r ( ft )}, h — T/n, 
i=i 

where only the variations under a given threshold function r(h) are taken into ac- 
count. The first main result of our paper is showing the consistency to ICt, as 
the number n of observations tends to infinity. Not equally spaced but synchronous 
observations are allowed for such result. The second group of results is given in 
presence of only FA jumps. If we dispose of non-synchronous data we still reach 
consistency by modifying our estimator in a similar way of [13] and [12]. When 
observations are evenly spaced, we prove a joint CLT delivering: 1. that ICt,u is 
asymptotically Gaussian and converes with speed yh, which extends results in [2] 
who estimated ICt in absence of jumps; 2. consistent and asymptotically Gaussian 
estimators of the regression coefficients /3s and of the correlation coefficient of the 
two continuous parts of processes 

In a further paper ([11]) we explore the speed of convergence of the estimator ICt,u 
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proposed here in the presence of, possibly correlated, infinite activity Levy jump 
processes and J^ 2 \ with dependence structure described by a Levy copula. 

The threshold criterion originated in [21] to separate the diffusion and the jump 
parts of a univariate parametric Poisson- Gaussian model. The criterion was shown 
to work even in nonparametric frameworks in [22], [23] and [15] . Potentially the 
threshold technique can be useful in each context where disentangling the quadratic 
variation of a signal has some importance to capture the contribution given by the 
diffusive component of the model and the one given by the jump component. 

The literature on non parametric inference for stochastic processes driven by dif- 
fusions plus jumps, based on discrete observations, is mainly devoted to univariate 
cases. As for bivariate processes Barndorff-Nielsen and Shephard Q4J) and Jacod 
and Torodov ([IH]) explore tests for the presence of co-jumps based on estimators 
constructed basically using cross multipower variations. 

We adopt the threshold method here since, at least in the finite activity case, it 
is a more effective way to identify (asymptotically) the intervals between consec- 
utive observations where jumps occurred. In fact already in the univariate case 
the threshold estimator of IV ^ is efficient (in the Cramer-Rao inequality lower 
bound sense), the asymptotic standard estimation error being a/2 x IQ^\ where 
IQ^ := Jq (a^) 4 dt (see [23]), while the multipower variation estimators are not 
efficient (the standard errors are all higher, see [3], [28] and the discussion in [23J, 
section 3.3, the infimum being a/2.609 x IQW). For the bipower covariation based 
estimator (BPC) of ICt, a CLT has been shown to hold only in absence of jumps 
([4j), in which case the standard error is a function of pt, <j\ ; , > which for instance 
equals 1.3 x J T (1 + p 2 t ){af ] fdt if (<xf } ) 2 = (a 4 (2) ) 2 , while here we show that a CLT 
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holds for the threshold estimator even in presence of (finite activity) jumps, the 
asymptotic standard error being J T (1 + p\ )((Tf) 2 (a^) 2 dt, for any p t , <jf \ and 
which is less than the error of the BPC at least when (cr^) 2 = (a^) 2 . The bipower 
covariation test of [I] has been discussed by j7], where the Authors show that, when 
dealing with large portfolios, it is necessary to use a different global cross-variation 
index to get reliable results. 

A CLT using multipowers for a bivariate process and in presence of jumps is given 
by [19j. More precisely, regarding the co-jumps, they consider the quantity B n := 
^^ =1 (A i X( 1 )) 2 (A i X( 2 )) 2 , which is an estimate of B : = £ s < T (AX s (1) ) 2 (AX s (2) ) 2 , not 
directly comparable with an estimate of the sum of the co-jumps 
we give here. Their goal is to give a test for the presence of co-jumps, so they 
concentrate on 

, W = E,(A^ (1) ) 2 (A^ (2) ) 2 
9n Ei(A i X«) 2 (A i X( 2 )) 2 ' 

the quotient of two cross-power variations of the bivariate X, computed for different 

lags kh and h: := X® - x[f_ 1)kh , A { X^ := X% ] - X { ( ?_ 1)h . They reach that 

0n ' ) — > 1 as n — > oo, on the space f2 < -^ where some co-jumps occur, and they prove 

a CLT for in restriction to We remark that to compute an estimator of 

the conditional asymptotic variance of (pi they in fact use the threshold technique 

when the volatilities are stochastic and are allowed to co-jump with the respective 

state variables. 

An outline of the paper is as follows. In section 2 we illustrate the framework; 
in section 3 we deal with the case where each component of X^ q ' has finite 
activity of jump. We show that ICt,u is asymptotically Gaussian, so that it is also 

6 



consistent. We find a joint CLT allowing to estimate the (3s and the correlation 
coefficient of the continuous parts of the two processes X^ q \ and we deal even with 
the case where we dispose of non-synchronous observations. In section 4 we deal with 
the more complex case where each can have an infinite activity semimartingale 
jump component We show that our estimator is still consistent. Since the given 
theory asserts that we can asymptotically identify the quantities of our interest, in 
section 5 we check on simulations that in fact the finite sample performance of ICt,u 
is good even for time step between the observations large enough (five minutes) to 
avoid considering microstructure effects on the data, at least for commonly used 
financial models with realistic choices of the parameters. Section 6 concludes and 
section 7 contains all the proofs and technical details. 

2 Framework and notation 

Consider a filtered probability space (Q, JF, (T t )te[o,T\, P) where X* 1 ) = (X^) te [ 0jT ] 
and = (JT t ^) t6 [ 0) T] are two real Ito semimartingales defined by 

Xl q) = f a<pds + f a^dW^ + 4 q \ t G [0, T], q = 1, 2 (1) 
Jo Jo 

where 

Al. = (W t (1) ) t6[0) T] and = (W t (2) ) te[0) r| are two correlated Wiener 

processes, with quadratic instantaneous covariation given by 

d < > t = p t dt, t G [0, T]; 

we can write 

= fHW® + y/T^ t W} 3 \ 
7 



where and are independent standard Brownian motions. 



A2. The diffusion stochastic coefficients = (a\ )te[o,T], 

= (a[ 9 ' )te\o,T], Q = 1,2, and p = (pt)te[o,T] are cadlag adapted 
processes. 

As for the jump components J( g \ in the next section we have FA jumps, i.e. 

k=l 

as specified with more detail below, where = (JVj )te[o,T] are counting processes 
with E[N^} < oo. 

More generally in section H] each is allowed to be any pure jump semimartingale 
with possibly IA. 

To begin with we assume to have equally spaced and synchronous observations. 
The consistency results under not equally spaced but synchronous observations are 
straightforward using Lemma 12.11 and Theorem 17.11 below. Generalization to not 
equally spaced and not synchronous observations are dealt with later. Let, for each 
n, 7c n = {0 = < < ■ ■ ■ < tn = T} be a partition of [0, T]. For simplicity 
let us write tj in place of Define h := tj — tj-i = ^, for every j = 1, ..,n and 
n — 1, 2, ... Note that h — ► if and only if n — > oo, so when computing the limits of 
our interest we indifferently indicate one of the two. 
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A3. We choose a deterministic function, h i— > r(h), satisfying the following 
properties 

/ x h log -r 

lim r(h) = 0, lim ,\ h = 0. 
h-*o h^o r(h) 

We denote r(h) by r^, and, for each q = 1,2, 

Dl g) = [ a( g) ds+ f ai q) dW^, 
Jo Jo 

the Brownian semimartingale part of X^ q \ 

As a consequence of the Paul Levy result about the modulus of continuity of the 
Brownian motion paths, we can control how quickly the increments of the diffusion 
part of each AjXW tend to zero. This is the key point to understand when AjX^ 
is likely to contain some jumps. More precisely, the Paul Levy law (see e.g. [20J, 
p. 114, Theorem 9.25) implies that 

a.s. lim sup — 3 < 1. (2) 

fc-°je{i,..,n} J2h\og\ 



However the stochastic integral a.W is a time changed Brownian motion ([27], The- 
orems 1.9 and 1.10), i.e. Aj (a.W) = Bjy t — Bjy t i , where B is a Brownian motion 
and IVt is the integrated variance J* a 2 s ds up to time t. Note that the increments of 



the drift part of X tend to zero more quickly than y 2h log | as h — > 0, so for 

we can reach a result similar to (j2J), as soon as the boundedness of the paths of a 

and a is guaranteed (which is the case when they are cadlag). In fact 

\J tj a s ds + p 3 a s dW s \ \ J '? j a s ds\ \ J '? j a s dW s \ 

sup — , < sup — j^=^=- + sup — — < 

i =1 »« \Mlog± >' v/2/ilogi i j2h\og\ 
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, HT l^-^J yj2A.IV log y^MMlog^ 

C{uj)aI r +sup sup , sup , 

V log ^ ' ^AjIVlog-sjjy i ^2M{u)h\o gj ^- h i ^/21ogi 

where C(o;) := sup se r T i |a s (a?)|, M(a;) := sup sg r T i |cr^(a;) | . By [20] (Theorem 9.25) 

and the monotonicity of the function xln ^ it follows that as h — > 0, the right hand 

side has a limsup which is bounded by ^/ M{u), thus for sufficiently small h, even 

in the case of not equally spaced observations, the following holds. 

Lemma 2.1. ([23]) Under A 2 we have that, given an arbitrary partition {t = 
0, ti, t n — T} of[0,T], then for sufficiently small h := sup - =1 n \tj — tj-i| we have 
a.s. 

I A I 

sup J < ^(w), g = 1, 2, 

i= 1 - ^2/1 log i 

where K q (uj) := a/ M(c<j) + 1 are finite random variables. 

Last result implies that if (AjX^) 2 > r^ and is, for small h, larger than 
2i^/tlog~ (as it is, under A3), then we have (AjX^) 2 > 2K 2 h\og^, and it is 
not likely that AjX^ coincides with the increment of a Brownian semimartingale, 
while it is likely that some jumps occurred within t,-] and made |AjX^| large. 

Application of Lemma 12.11 gives us the main tool for the construction of our 
estimators in the next section. 

Notation. 

• For any semimartingale Z, AZ S = Z s — Z s _ denotes the size of the jump of Z 
at time s, while AjZ = Z tj — Z t] _ 1 denotes the increment of process Z in the time 
interval tj] 

• IC t = J * PsCTs 1 ' ] (Ts ] ds denotes the integrated covariation up to time t, 

IC t>n = J2j=l..n: t s <t A i- X ' (1)l {(A J -XW)a<r h } A i- X ' (2)l {(A j XW)!»<r k }' h = T / n ' IS its thresh- 
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old estimator 

• IV t ^ = f*(o~i 9 ^) 2 ds denotes the integrated variance of process q=l,2, up 
to time t and IV^ n = n - t j <t(^-j^-^) 2 ^{{A j x('i'>) 2 <r h } is its threshold estimator 

• sometimes AjX^l^ AjX ( q )^2 <rh y is indicated briefly with AjXi 9 ^ 

st 

• sometimes we write Plim to indicate the limit in probability. — > indicates 
stable convergence in law of processes. See [18] . ch. 8, sec. 5c, for the definition 
and properties of stable convergence in law, and [16] for further statement of useful 
properties 



3 Finite activity jumps: consistency and central 
limit theorem 

In this section we assume that is any FA jump process: for each q = 1,2, 

4 q) = / tW^Ev 

Jo T~, k 



k=l 



where = (N^ )te[o,T] is a counting process with E[N^] < oo, {r^ , k = 

1, jVy } denote the instants of jump of J*- 9 - 1 and 7 ( 9 ) denote the sizes AJ t (9) of the 

r fc 

jumps occurred at rj. . Denote 

7 (9) = min 

~ k=l,...,N^ q) k 

A4. Assume E[N$] < 00 and Ph w = 0) = 0, V jfe = 1, q = 
1,2. 



Remark 3.1. Condition A4 implies that a.s. y( q > > 0. 
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Example 3.2. If are FA Levy processes, then they are of compound Pois- 
son type ([9], Proposition 3.3, section 3.2): are simple Poisson processes with 
constant intensities and for each q the random variables 7 (,) are i.i.d., for 

T k 

are independent on N™' and satisfy condition A4. 

We remark that the consistency and CLT we reach in this section are valid 
in presence of general finite activity jump processes, in that we do not need any 
assumptions on the law of the jump sizes, or of the counting processes nor 
any assumption of independence. We do not even need that are FA jumping 
semimartingales, we only need that A4 holds, which is true if are (FA jumping) 
semimart ingales . 

Now we construct our threshold estimators. 

Definition 3.3. We define for r,l 6 IN 

j:tj<t 

wM(xw,xv>) t = h- 1 n A ;+^ (1) ri A ;+^ (2) - 

and their analogous threshold versions 

~W (x (i U ( 2))< = h i-r±l (A,A«)^l {(AjX(1))2 < r4 (A,A( 2 )) i l {(AjX(2))2 < r4 , 

j:tj<t 

1 1 

w^ix^x^h- 1 ^2 nv^ (i \A,,x( 1 »p< ril n^ i(2)i f(A,,^F<,}' 

vlf(XW,XW) T and w^ n \X^ l \X^) T are used in [2J to estimate IC T in the case 
where X^ are diffusion processes, vjw (X^ , X^)x and (X^\ X^)t are mod- 
ified versions for the case of jump- diffusion processes: by Theorem 1 7 . 1 1 1 hey exclude 
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from the sums the terms containing jumps. Note that v^i(X^\ X^) t = IC t , n , for 
all t E [0,T]. 

In view of the practical application of our estimator we are now interested in the 
speed of convergence of ICt,u- We in fact reach even more. The first main result of 
this section is a joint central limit theorem for the threshold estimators 

1 fv {1) IC X 



V 



IC IV 



(2) 



which implies that in presence of finite activity jumps /CV,n converges to IC at speed 
Vh, h = T/n, and it allows to give estimators of standard dependence measures 
between the diffusion parts of our processes X^ 9 \ such as the realized diffusion 
regression coefficients up to time t 



(1,2) .. 



ic t 



IV, 



(2)' 



(2,1) 



ic t 



IV 



(1) 



and the realized diffusion correlation 



Pi 



(1,2) .. 



ic t 



Theorem 3.4 (Joint CLT, FA jumps). Under assumptions from Al to A4, with 
h = T/n, we have, as h — > 0, 



( IV™ - IV (1) IC n -IC \ 



\ 



I 



IC n - IC IV® - IV {2) 
where Z is the 2x2 process with components 



St J_ 



2Z\i Z\2 + ^21 

Z\2 + Zi\ 2^22 J 



Z lu := / (o-™) 2 dB lls 

'0 
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Zi2,t ■= Jq p s Os<?sdB lu + /J a/1 - pi a (2) dB X2s 
Z 2 i,t ■= J t Ps ^ 1) ai 2) dB lls + y/T=fi ai 1) ai 2) dB 21s 

Z 2 2,t := j\ 2 MPfdB ns + £p.y/T=fi (of ) 2 [dB Us +dB 2U )+j\l-p 2 s )(a^) 2 dB 21s 

and B is a 2 x 2- dimensional standard Brownian motion independent on the filtered 
probability space te[o,T]> P) where our model is defined. 

Note that the result for IV^ is consistent with E3L since 



Var(V2Z qq , T ) = 2 / 
Jo 



T pT pT 

4/_(o)\4j„ , a / J/1 „2w M\4 , . / /, „2\2/_fo)\4 



= 2 r {o-^yds. 







Corollary 3.5 (Consistency, FA jumps). Under Al to A4, as n — > oo ; /or a// 

t G [0,T] 



iCt.n ICt, 



if a.s IV t 3 ^ i/ien 



A { ^:=^-# J) , (*,;;) = (1,2), (2,1) 

41,2) _P (1,2) 

' / : (i) : (2) ^ ' 

Corollary 3.6 (Speed of convergence of /5s and p, FA jumps). If a.s. IV ^ ^ for 

all t G [0, T] we have, for = (1, 2) or (2, 1), 

7,-1/2 ( _ ^ Z 12 + ^21 ,Hy 7 IC 
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// o.3. IVfhv? ^ for all t e [0,T] we have 

IC 



h- 1/2 P^ 2) 



St 



Z\2 + Z 2 l Z22IC Z\\IC 



The following proposition allows us to give a CLT for the standardized version of 
the estimation error ICt,u—ICt- Note that the asymptotic variance of /i~ 1//2 (/CV,n — 
IC T ), by Theorem^! is given by (Var(Z 12jT +Z 21 , T ))/2 = /^(l+p 2 )^) V t (2) ) 2 dt. 

Proposition 3.7 (Estimate of the standard error for IC t>n , FA jumps). Under 
assumptions Al to A4 we have, for all t e [0, T], 

vi n 2 \x^\X^) t -w^(X^\X^) t -A Al + ^)(a«) 2 (af) 2 ^. 



Jo 

We now are ready to present the central limit theorem for the standardized estima- 
tion error. 

Corollary 3.8 (CLT for the standardized version of ICt, n — ICt, FA jumps). Under 
Al to A4, if a.s. + p 2 )(ai 1} ) 2 (af ] ) 2 ds ^ we have 



^h\j4i ,X (2) )(-» (n) {XV , A( 2 ) j t 
where M denotes a standard Gaussian random variable. 



/V, 



Remark 3.9 (Estimate of the co-jumps). By Corollary 13.51 clearly we have an 
estimate of the sum of the co-jumps up to T simply subtracting ICt,h from the 
quadratic covariation estimator: 

J^AjX^AjX^-ICt^^ M 1} AJi 2) , 

3=1 0<s<T 
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as n — > oo. Analogously we can obtain an estimator of the sum of the co-jumps up 
to each time t e [0,T]. 

An estimate of each AJ^ AJs 2 \ with s 6 [0, T], is obtained using 

AjXWAjX® - A J A( 1 )l {(AjX(1))2 < ra A J X( 2 )l {(AjX(2))2 < rh} , (4) 

with j such that s E]tj_i,tj]. Alternatively, as we consider one single term, and not 
the sum of n terms, even 

A i X {1) l {(AjX (i) )2>r;i} A i A (2) l {(AjX (2) ) 2 >rh} (5) 

or 

AjX^AjX® (6) 

estimate the co-jump AJg AJs 2 \ with s &]tj-i,tj], since | AjX^l^ A . X ( q )^2 <rh yAjX^ 
< 2 v ^sup sG[0jT] \X(% q = 1,2, £ = 3 - q, and lA^Wll^w^^jlA^^lx 
!{(AjX( 2 )) 2 <r h } ^ r ft tend to zero in probability as h — > 0, by the pathwise bounded- 
ness of each X^> on [0, T\. However, as we show in section 5, estimator (jHJ) has the 
best finite sample properties in the simulations of Model 1 having FA jumps. □ 

Remark 3.10. Finite sample performance and micro-structure noises. Our 

theoretic results allow to estimate ICt and the co-jumps asymptotically for h — > 0, 
while in practice for very small values of h financial time series are affected by 
microstructure noises which introduce a bias which is larger as h is smaller. In 
section 5 we implement our estimators of the integrated covariance and of the co- 
jumps on simulations of realistic financial time series and we find that they have good 
performance already with temporal mesh h corresponding to five minutes, a time 
lag at which prices are not usually affected by microstructure noises ([5]). However 
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we remark that when iid micro-structure noises contaminate the observations of each 
asset price X^ q \ the threshold estimator rules out even the noises, similarly as it 
rules out the contribution of the jumps ([21]). 

Remark 3.11. Asynchronous observations. It is known that the problem of 
the estimation of the covariation among two assets undergoes the so called Epps 
effect, i.e. in the empirical applications the estimator tends to zero as the step of 
observation h tends to zero. The asynchronicity among the observations of 
and X^ is considered one of the possible causes ([26]; [5], section 2.10.2). In fact 
some Authors have tackled the problem of reaching a consistent estimator of the 
covariation even when data are asynchronous and h — > 0, under the assumption of 
Brownian semimartingale models (in [13] the estimator is introduced, however we 
refer to P2] where the observation times are allowed to be dependent on X^ q '). 
At the time scale of five minutes the Epps effect probably does not affect our estimate 
of ICt- However even in presence of this microstructure-type noise (for smaller h) 
it is possible to make our estimator correctly converge to the integrated covariation, 
as detailed below. 

Assume we dispose of two records {D {l f n) , D^( n) , ...D^( n) }, {D™ n) ,D( n) , ...D^ n) }, 
of observations of two Brownian semimartingales D^ 1 ' and D^ 2 \ with the two stochas- 
tic partitions = Tq < rj^ < ...r and = < < ■■■v^l) spanned on 
[0,T]. For simplicity let us write z/, and tj in place of vf 1 ^ and r,-"\ The idea 
of Hayashi and Kusuoka is to select only some of the cross variations (D^ — 
Dr}_ 1 ){Dui — -Di>i_i), in order to estimate the covariation, and precisely the ones 
for which there is an intersection between the time intervals ]tj_i, Tj] and ]vi-i, vi]. 
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We show here that using their result ([12J, Corollary 2.2) we in fact reach the same 
kind of consistency in the case of asynchronous observations even in presence of 
finite activity jumps. The idea is very simple: we first eliminate the jumps, us- 
ing threshold technique, and then apply the Hayashi and Kusuoka estimator to the 
estimated continuous components D^ q ' . Recall that 

(x^ - xl% = - ^ (9) )i {( 4,)_ x ( 9 ) )2 < rfc} , 

for any two time instants a and b, h := sup J=1 m ( n ) (tj — Tj_i) V sup i=1 >>fc („) — fi-i). 
We in fact have the following 

Theorem 3.12 (Asynchronous observations). Let Al to A4 hold, = r < t% < 
... < T m {n), = Vq < V\ < ... < u k ( n ) be two sequences of stopping times such that 
r m (n) t T, z/ fc( „) | T a.s., asn^oo then 

p 

as h := 8ixp j=Um w(Tj - 7y_i) V sup i=1 ..*(») (fj - -> 0. 

4 Infinite activity jumps: consistency 

In this section we allow the jump components of processes to have infinite 
activity, so we are here in the case where are general Ito semimartingales. 
Any unidimensional Ito semimartingale has a representation as in ([1]) with each 
decomposed as 

= JoJ h (< H .,t, x) \>i^ q Ku,t,x)^\u,dx,dt), (7) 
•4 ? V) = j?J] 7 ( 4 ) (ta , M )[<i7 (ff) ( w )*> a: )A (,) ( w > da: » tfe )> 
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where fj,® is the Poisson random measure of the jumps of j( q \ ^ q \uo,dx,ds) = 
dx, ds) — dx, ds) is its compensated measure, j^ q \uj, dx, ds) = dx x ds, 

the coefficients 

a (<z) ;(T (?) )T (9) are 

predictable and JlA (j^) 2 (u,t,x)dx is a.s. finite 

(see H6], pp.3,4; [IS], (2.11)). 

Conditions A2 and A4' below guarantee local boundedness properties of such co- 
efficients. 

A4'. JlA {j ( - q) ) 2 {u,t,x)dx is locally bounded. 

For each q = 1,2, j[ 9 ^ is a finite activity jump process of type j[f = Ylk=i 7 
as in section [31 where E[N^] < oo is equivalent to £7[ f («)i >:L dx] < oo and now 
the sizes I7 (,)| are all larger than 1; on the contrary accounts for the infinite 
activity jumps of J^ q \ since generally Ji^g)^ dx = +00. is a compensated sum 
of jumps, where each jump is bounded in absolute value by 1. Therefore, for each 
q = 1,2, j[ 9 ^ accounts for the "large" and rare jumps of X^ q \ while J% accounts 
for the frequent and small jumps. 

Example 4.1. If one of the two processes is a pure jump Levy process, it is 
always possible to decompose it as in ([7]) with 7™(w, t, x) = x but v^ q \u), dx, ds) = 
u^ q '(dx) x ds, where is the Levy measure of J (9) and is a deterministic a-finite 
measure such that J u 1 A x 2 ^ q \dx) < 00 but generally such that f^ <x v^ q '{dx) = 
+00. 

We prove that ICt, n is still a consistent estimator of ICt, for all t G [0, T]. For 
ease of notation we only consider IC up to time T and evenly spaced synchronous 
observations. Not evenly spaced but synchronous observations (with h = sup^- \tj — 
tj-i\) and arbitrary t G [0, T] are straightforward. As a consequence the same 
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estimators of the co-jumps, presented in the previous section, are consistent even in 
the present framework. 

As for the speed of convergence of ICt,u, in the presence of infinite activity jump 
components, things are more complicated in that such a speed is determined both 
by the dependence structure between J2 an d by the amount of jump activity 
of each J% . In [H] we consider two Levy infinite activity jump components j% and 

~(2) 

Jr, with a dependence structure described by a Levy copula. We find that, when 
and do depend, the speed is still \fh only when the activity of jump of at 
least one process is moderate (Blumenthal-Getoor index smaller than 1), otherwise 
the speed is less than \fh. 

We now state the main result in presence of infinite activity jumps. 

Theorem 4.2 (Consistency in presence of IA jumps, synchronous observations). 
Let (X' t ;)te[o,T] an d {X\ )te[o,T] be two processes of the form (Qp. Assume Al, A2, 
A3 and A4'. Then 

p 

ICt,ti — > ICt, 

as n — > oo. 

Remark 4.3 (Estimate of the co-jumps). Even in this framework of infinite activ- 
ity jumps, as a consequence of Theorem I4.2[ the sum of the co-jumps up to T is 
estimated by 

n n 

{(A J X( 1 )) 2 <r h } A i^ (2)l {(A J X( 2 )) 2 <r h }- 

i=i j=i 
Estimates of a single co-jump AJg AJg 2 \ with s £}tj-i,tj], exactly as in section 
[31 are given by (jlj), ([S]) or Simulations in section 15.31 show that for Model 
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2, with IA jumps, in fact estimator (JSJ) is a little bit more biased than ([6]) but 
is still acceptable. Note that since in sec. 5 each has infinite activity and 
j( 2 ) = pjJ^ + a/1 — p 2 jJ^\ each }tj-i,tj] contains an infinite number of co-jump 
instants. □ 



5 Implementation 

5.1 Choice of the threshold 

Our estimators depend on the threshold function r^. In this section we check on 
simulations how the results are sensitive to the choice of in a given class. This is 
only an informal and necessarily limited investigation. Formal study of methods for 
optimal threshold selection in a given model is object of further research. 
In principle there are many functions satisfying conditions A3. However on 
simulations we find that the choice of within the family of powers of h, = ch^, 
with c a constant and (3 a power in ]0, 1[, seems to be sufficiently good. 
We simulate two kind of models: Model 1, proposed in [13], where each has 
stochastic volatility and a FA Compound Poisson jump part and Model 2, proposed 
in [8], where each has constant volatility and IA jumps, as described in Table 
A. For Model 1 the parameters of the univariate X^ are taken from [TJ]. A path of 
each a varies most between 0.013 and 0.019 in a day. For Model 2 the parameters 
of the univariate XS^ are taken from Table 2 of [Sj for GE and HWP stocks. Note 
that the parameter Y is not significantly different from zero for the two considered 
stocks, so that the CGMY process can be reduced to the VG process. 

The VG process is characterized by three parameters k, 9 and <;. It is obtained 
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by evaluating a Brownian motion with drift, 9t + <;B t , at a random time G t given by 
a gamma process, a Levy process whose lag h increments G t +h — G t are distributed 
as Gamma r.v.s with mean h and variance Jik. It turns out that the VG process 
is pure jump and has infinite, but moderate, activity (it is a process with finite 
variation). 

To effectively introduce non zero co-jumps, in each model the jump component 
j( 2 ) of is correlated with of in the following way: we generate and 
an independent with parameters as in Table A, then = pjjW + y/T^ptj®. 
The simulation of the model paths has been made using the Euler scheme with 
increments of 1 second, then we have taken the five minutes synchronous returns 
and constructed our daily threshold estimator ICt,u- We simulated 3000 bivariate 



For each model we implement the estimator of IC as varies. Figures 1-2 show 
how the mean relative bias in percentage form 



varies as (5 varies in ]0, 1[ for c = 0.1, ....,5.6 with step 0.5 in Models 1 and 2, for h 
fixed equal to five minutes (n = 84 observations per day, time unit of measure T—l 
day, h = 1/84). It is evident that the choice c = 0.1 is the best one since in presence 



the case of IA jumps the bias is much larger but c = 0.1 allows to reach, for high /3, 
the lowest possible error. Figures 3-4 show the empirical densities and the QQ-plots 
of the normalized bias 



paths. 



100 



(IC Ttn - IC) 
IC 



of FA jumps (A^ = 0.118) it allows to decrease the bias as (3 increases. In fact in 



ICt,u — 



IC T 



(8) 
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when r h varies as before, for fixed h equal to five minutes, for Model 1 with A^ = 
0.014. The same plots for Model 1 with A (<?) = 0.118 and Model 2 are shown in 
Figures 5-6 and 7-8 respectively. We conclude that the best choice is r h = 0.1/i°". 
As a further check in Figures 9-10 we made the same plots for Model 1 with \® = 
and we found that the choice of r h gives good results as well. 

5.2 Estimates of IC and Eo<i<T A J^^J? ] 

on simulations 

We report here the performance of the estimators of ICt and of the sum J2o<t<T 
of the co-jumps up to time T, where the threshold is the one selected in the previ- 
ous subsection. T is kept fixed to one day, h equals five minutes. Figures 11-12-13 
show the histograms of 100 ( ' ICt i " ( 7 IC ^ to check the efficiency of IC T , n f° r Model 1, 
A (<?) = 0.014, Model 1, A (<?) = 0.118 and Model 2 respectively. The relative summary 
statistics are shown in Tables 2 and 3. ICr,n has an acceptable performance in 
Model 1 and it is biased in Model 2 but note that the estimation errors for IC T 
and J Q (1 + Pt)(<Jt ){pt ) 2 dt compensate and give good empirical densities of the 
normalized bias in Figure 7 for c = 0.1 and (3 = 0.99 and Table 1. Figures 14-15-16 
show the histograms of the following relative bias in percentage form 
1 nn (E? = i A,XW A,X( 2 ) - IC T , n ) - Eo<t<T AX f (1) Axf } 

— (I) — 

for the sum of the co-jumps in Model 1, A^ = 0.014, Model 1, A^ = 0.118 and 
Model 2. Note that since in Model 2 each is a pure jump process with IA 
in fact each movement of is a jump, and each time that jumps even 

J(2) 

does by the way we correlated them, therefore the best we can do to reach the 
true Eo<t<T Axf 1 ' 1 AX^ 2 -* is to take the sum of the cross-products of the one-second 
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differences of processes J^ q \ Tables 4 and 5 show the relative summary statistics. 
The performance of the estimator of the sum of co-jumps is very good under Model 
1 and a bit worse under Model 2. Under model 2 the estimate of ^ 0<i<T AJ^AJ^ 
is much better than the one of ICt- 

5.3 Estimate of the single co-jumps 

Using the threshold function selected in subsection 15.11 both for Model 1 and for 
Model 2 we implement (jSJ) and (jSJ) to estimate each single co-jump, in order to 
check which is the most informative estimator for the single co-jumps. We consider 
1 day time horizon and h equal to five minutes. Figures 17-18 show the histograms 

JJ -AJ (1) AJ (2) 

of the 3000 values of 100 — - — — jrf— f° r each estimator for both Model 1 (Ai = 
A2 = 0.118) and Model 2, where we define by JJ t (joint jumps) the estimate of 
A Jt AJt . Tables 6-7 report the relative summary statistics. We conclude that 
the most informative estimate of AJ t A J t is ((SJ) for Model 1 and t^yfj for Model 2, 
however, we find that (El) for Model 2 is still well acceptable. We find that, anyway 
all three estimators show a good performance, since the mean percentage estimation 
error in the worst case (estimator (4) Model 2) is 1% with low standard deviation. 

6 Conclusions 

In this paper we introduce a new estimator of the diffusion part IC and of the 
co-jumps in the quadratic covariation of two semimartingales X^ q \ To capture the 
separate contributions to the quadratic covariation has important applications in 
finance (forecasting, option pricing, risk and credit risk management). 
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The estimator ICr, n is constructed using a threshold criterion introduced in [21], 
and consists in summing properly selected cross products of increments of the two 
processes. Our estimator is consistent, and when the two jump parts have only finite 
activity a joint CLT for ICt,u and the estimators IV? of the integrated variances 
is proved and delivers the following important consequences. 

1. ICt,ti is also asymptotically Gaussian with speed of convergence \/h. A central 
limit theorem in presence of infinite activity jump parts is studied in a further pa- 
per ([H]) where we find that the speed of convergence of /CV, n is determined both 
by the dependence structure between the two processes X^ q ' and by the amount of 
jump activity of each J™. 

2. Consistent estimators both of the sum of the co-jumps occurred within [0, T] and 
of each single co-jump are obtained. 

3. We construct asymptotically Gaussian estimators of the regression coefficients j3s 
and of the correlation coefficient between the two processes X^ g \ 

Further we find that in presence of FA jumps a slight modification of ICV, n is con- 
sistent even when only non- synchronous observations are available. 
We assess the choice of the threshold and check the performance of our estimators on 
two different kind of simulated models which are common in the financial literature. 
Model 1 has components with stochastic volatilities and FA jumps, while Model 2 
has components with constant volatilities and IA jumps. We find that even with 
five minutes observations the performances of the estimators of ^2 0<t<T A AJ^ 
and of the single co-jumps are satisfactory. ICt,u is satisfactory in Model 1 while is 
biased in Model 2 but the corresponding normalized bias has still Gaussian behavior. 
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7 Appendix 

The following theorem is the key result, in the finite jump activity case, validating 
the idea that if (A^X^) 2 is larger than then some jumps occurred in ]tj_ 1 ,tj] 
(and vice- versa). It is stated in the general case of not equally spaced observations. 

Theorem 7.1. (JEBj, FA jumps) Under the assumptions from Al to A4, given an 
arbitrary partition {t = 0,ti,...,t n = T} of [0,T], then for sufficiently small, but 
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strictly positive, h := sup =1 n \tj — (depending on u) we have a.s. 



l{(AjX(«)) 2 <r fe } — 1{AjJV(«)=o}j J — 1 ; 2, n, g — 1,2. 



□ 



Proof of Theorem \3-4\ [Joint CLTJ By Theorem 17.11 we have, for all 

te [o,T], 



t,-<t 



,-iVC 1 ) =0^3^(2) =0} 



/a 



.A^l, 



ti<t 



.ivw^o} Aj£) (2) + /i 1/2 Aj-D (1) l{ Aii v(i)^ 0} ^i^ 'i {AjiV (2)^ 0} . 



-A-V^A^Wl^ 
Each one of last three sums tends a.s. to zero as h — > 0, since it contains at least 



one I{AjNM^o} an d fc> r any q = 1,2 we have 
^^A^WA^V 



Plim 

71— >00 



< P\imK 1 (cu)K 2 (uj)Vh\og]-N^ ) = 0. 



Moreover, analogously as in [23], for each g = 1, 2 we reach that 



h- 1/2 £( A ^ (9) )* - Ivi t 9) = h ~ 1/2 5>^) 2 - iv\ 



(?) 



h- 1 * j>^ (9) ) 2 -^ 9) -h-^J2(^ D(q) )%^^ 
\tj<t j tj<t 

where the last term tends a.s. to zero as h — » 0. Therefore we have that 



1 iv { y 



IV {1) ic T>n - IC 
(2) IV^ 



IC n - IC IV n 



J 



(9) 
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has the same limit in distribution as 

£,,< t (A^M) 2 _ jy« AjDMAjD® - IC > 



Note that 

^ AjDMAjD® -IC t = J2 (A^WAj-DW - A, < D^-\D^ > ) (10) 

and, along the lines of |5| (proof of Theorem 1, sec. 3.1), using Ito formula we know 
that 

d(D^D^) = D^dD® + I^W 1 ) + d < £>« D (2) >, 



so 



Therefore 



AjD^AjD^ = Aj(D^D^) - D^AjD® - D^AjD^ 

= f D^dDf^ r Df]dD^ + A ] <D^\D^>-D^A ] D^-D^_ i A 3 D^\ 
so that ( fTUl) equals 

v t<t ' ./o v 7 



i j(n) 



where 



^(n) 



* n. ' * n ' 



(2) 
s 
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As special cases, for each q = 1, 2 



^(A^) 2 - IV { t g) = ((A^ (g) ) 2 - A? < ^ ^ > ) 



t n <t 



2 f (D® -i2 D ^h^ m {s))dD? = 245, 
By Theorem 5.5 in ^T7\ we have that 

h -l/2 A (n) ^ Z 

with Z as in (J3]). It follows that, as n — > 00, (JSJ) converges stably in law to 



1 

7! 



^ 2Zj 1 Z 12 + Z21 ^ 



V 



/ 



□ 



Z\2 + ^21 2^ 2 2 

Proof of Corollary \37U\ [Speed of convergence of (3s and p, FA jumps] 

For all t E [0, T] we have 



h- i/2 m> 



therefore 



IC \ 



IV U) J V2IV ij) (IV®) 2 ' 
As for pn ,2 \ note preliminarily that Theorem 13 . 41 implies that h~ 

-1/2 f ^iv® - J IV 



converges stably, since, t by t, 



( jyU) _ fy® 



si 
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As a consequence 



ic \ 
VrWrv^) 



rti) 



V2 IV® 
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h _ 1/2 IC n IC _ + h - V 2 jn 



Jiv {1) fv {2) \Jiv {1) iv {2) s/iWivv 

V n n \y n 71 

The first term converges stably to Zi 2 +z 2 ^ = ^jig fa e second term equals 

a/2 IVWIV^ 



h-^IC I l x/IV® l 



'ir [ -1C I I 1 \ h-V 2 IC 1 1 \ !Vr ] 



st Z22IC Z\\IC 



V2 iv^iiv^y/ 2 V2 iv {2 \iv w f/ 2 

□ 

Proof of Proposition \3.7\ [Estimate of the standard error for IC n , 
FA jumps] For t = T it is sufficient to show that as n — > 00 

v^\X^) T -A / T (2 Pt 2 + l)K«) 2 (^) 2 dt, 

./o 

and 







For t < T the proof is analogous with ^™ =1 replaced by ^2j. t <t • By Theorem 17.11 
we can write 



Phm vf}{X^\X^) T = Plim /.- 1 ^(A,D«) 2 1 {A37V(1)=0} (A^( 2 )) 2 1 {A 



Plim i^p^D^-Plim /T 1 V (A, j DW) 2 (A, j D( 2 )) 2 1 |a . 



n— >oo ' ' J 
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-Plim fr 1 J](A, j D( 1 )) 2 (A jj D( 2 )) 2 1 {AjJV( ^ 0} 
i=i 

n 

+5™ ^ E( A ^ (1) ) 2l {A^(^0}(A^ (2) ) 2 l{A^( 2 )^0}- 

By Theorem 2.1 in [6], 

Plim v^\D^) T = / T (2p 2 + l)(a; i) ) 2 K (2) ) 2 dt, 
whereas the other terms are all zero. In fact for any q — 1, 2 



Plim Jr 1 V(A, J D( 1 )) 2 (A, J D( 2 )) 2 1 {A ^ (9) ^ 0} < Plim K*(u)K*(u>)h( log = 0. 

(11) 

Now we deal with (X^\ X^)t- Analogously as before 

Plim w {n \X {1 \X {2) ) T 



n-l 



Plim h' 1 ^ 



n A i+i DW(i - i {(W ( 1¥0} ) n a ^ (2) (i - 1{(a j+i m^o } ; 

_i=0 i=0 

which coincides with the sum of Plim w^ n \D^\ D™)t with a finite number of terms 
which are shown to be negligible. By Theorem 2.1 in j6], Plim tfl^p' 1 ', D^)t = 
Jo Pt( <T t 1 ^) 2 ( a t 2 ^) 2 dt^ while the other terms are given by the product of Yl] =0 Aj +i D^ 
[]| =0 Aj +i £)( 2 ' with at least one of the indicators 1{a j+s at(<j)^o}' l° r an s G {0,1}. 
Therefore the limit in probability of each such term is zero as in (11 II) . □ 



Proof of Corollary \3.8\ [CLT for the standardized version of IC t , n 
ICt, FA jumps]. By Theorem 13.41 we have 



h- 1 ' 2 (fC n - IC 



St 



1 

7! 



a® a®. 



I-P 



dB 1 / + dB s 
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The variance of the last term at time t is J *(l + p 2 s )(a^) 2 (ai 2 ^) 2 ds. By Proposition 
13.71 we then obtain that 



where M is a standard Gaussian r.v.. □ 



Proof of Theorem \3.12\ [Asynchronous observations] Note that we can 
assume that a and a are bounded on [0, T] (|15j). so that the Brownian semimartin- 
gale parts of belong to L 8 . Using Theorem 17.11 in the not evenly-spaced 
observations case (|23j) with h := sup,- =1 .. m («)(Tj — t,_i) V sup i=1 (fc'i — Vi-i), a.s. 
for sufficiently small h we can write 

j=l..mW, j=i..fc(«) 

E (^-^J^^-iv^M^^ 

j=l..m( n ), j=l..fc( n ) 

E - ^)(^? - 

j=l..mW, i=l..fc(") 



j=l..m("), i=l..fc00 



^'-^^o} + 1 {iV«-iV« 1 ^0}+ 



1 



The first sum of the r.h.s. tends to ICt in probability by Corollary 2.2 in (|12j). 
with / = g = 1, while each sum in the second term is dominated in absolute value, 
for a suitable q, by 

bupI^-^JbupIDW-DWjjv*), 

which tends a.s. to zero as h — > 0, by Lemma [2.11 □ 
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The following facts are used within the proof of Theorem 14.21 
Without loss of generality (as in [15], Lemma 4.6) we can assume that 

A5. f x&R 1 A (7 (9) ) 2 (cu, t, x)dx is bounded. 

Lemma 7.2. For each q = 1,2 we have the following. 

1. If processes a and a are cddldg then, under A3, a. s., for small h, l{(A j D(<?)) 2 >r,j = 
0, uniformly in j . 

2. Under A5 we have that, for each j = l,..,n, E^AjJ^) 2 } < Kh, for a positive 
constant K . 

Proof. Part 1. is a consequence of Lemma [2. 1[ 
Part 2. 

EUAjJ^) 2 } = E\ f 3 [ (^f^^dx.ds)} = E{ [ (yto) 2 dxds] : 

since, by assumption A5, Jj 7 (, 7 )| <1 (7 ( ' <5, ' ) ) 2 ^ is bounded, the last term is dominated 
by Kh for some positive constant K. □ 

The following lemma generalizes analogous results given in [23] from the frame- 
work of Levy jumps to the one of Ito semimartingale jumps. 

Lemma 7.3. The following facts hold. 

1. Let us consider any sequence 7r n of partitions {0, t\, ..,t n = T} of [0, T], n G IN, 
such that maxj = i,. n \tj — tj-i\ — ► as n — > oo. For each q = 1,2, as long as 
is a semimartingale, we can find a subsequence for which a.s., for 
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any S > there exists a sufficiently large k such that for all j = 1, ..,rik on 



{(AjJ^) 2 < 4r h J we have 



(Ajg) 2 < Ar hk + S, Wse^t, 



2. Under A3 and Ab , for each q = 1,2, we have J2" =1 P{ AjN® ^ 0, (A, J 2 (g) ) 2 > 
Ar h } -> as /i -> 0. 

Proof. Statement 1 is a consequence of the fact that ([25], Theorem 25.1) 
there is a subsequence such that, defined h k = T/n k , ^^^(AjJ^) 2 tends to 
^ sg r 0t i(AJ2 9 j) 2 a.s. uniformly w.r.t. t G [0, T], as — >■ oo, where [x] denotes the 
integer part of x. Since a.s. 



sup \{AjJ\ 

j=l..n k 



(g)^2 



E ( Aj ; 

se]tj-_i,t 3 -] 



(<?)\2| 
2,sJ I 



{[tj/h k ] lti-i/h k ] 
[ J2 (a,j 2 ) 2 - E (Ajg) 2 ] - [ J2 (^ 2 ) 2 - E ( Aj S) 2 ] 
<=1 s6[0,%] *=1 »6[0,tj_i] 

[t/fcfc] 

< 2 sup [ ^ (A, J 2 W ) 2 - E ( A ^i) 2 ] - 0- 

we in fact have that a.s. for all j = l,..,Uk each squared increment (AjJ^) 2 is 
uniformly, on j, arbitrarily close to 5^ a e]t,--i *j](A^s) 2 - More precisely, a.s. for all 
o" > we can find a sufficiently large such that 



sup 

j=l..n k 



(A,J : 



(ah 2 



E ( Aj ; 



(<?h2 
2,s J 



so, for all j such that (Aj J^) 2 < 4r^ we have 



(Ajg) 2 < sup 

i=l..n* 



E ( Aj ; 



'(?)\2 
2,s J 



s£]tj-i,tj 



+ (A 3 J^) 2 <4r h + 6. 
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In particular for any s &]tj-i,tj] with j such that (AjJr, ) 2 < Arh, each squared 
jump size (AJ^j) 2 is bounded by 4r^ + 5. 

Statement 2. The predictable compensator of = J2 s <t^\Aj {q) \>i anc ^ ^ e P re " 
dictable quadratic variation of are of the form = J Q * As ds and A'} q ^ = 
Jq \'s ds respectively Assumption A5 guarantees that both Jj 7 ( 9 )| >:l 1 dx and 
Jj 7 ( 9) | <1 (7^) 2 dx are bounded, and therefore that A^ and A'^ are bounded pro- 
cesses. Using exactly the same argument as in [Tj, eq. (60), with 5 = 1 and ( = 3y/¥h 
and replacing M{8) = J j^^x^dx.dt) with our Jj , we conclude that 

^ P{AjN^ + 0, (Aj J 2 (9) ) 2 > 4r fc } = 0(n/i -). □ 

For any 5 > denote by Z^' 6 the following pure jump plus drift semimartingales 
having only jumps bounded in absolute value by y^r~h k + 5, g = 1, 2: 

:= f f j {q) ^ q \dx,ds) - [ [ 7 (<z) dxdt, t>0. 

By Lemma 17.31 we have that for any 5 > 0, for sufficiently large k the indexes j 
for which (AjJ^) 2 < 4r/ lfc are such that the increment (AjJ^) 2 coincides with the 
increment (AjZ^' 5 ) 2 of Z^' S , since (Aj J^) 2 does not contain jumps larger than 
\/±r hk + 8. 

Lemma 7.4. For each q = 1, 2 

n 

Plimy(A 7 ^' ?) ) 2 l / . A ?(9) , 2 ^ =0 

Proof Consider the sequence of partitions 7r n = {0, T/n, 2T/n, .., T}. Take 
any subsequence 7i nr By Lemma |7T3| point 1, a.s. there exists a sub-subsequence 
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ri£ k such that for any 5 > and k sufficiently large then for all j = l,..,ni k on 
(AjJ^) 2 < Ar hl we have (A.-J^) 2 = (AjZ^' 5 ) 2 . Denote 



n 



Therefore 



3=1 



0<Phm^) <Phm^(A,4f ) 2 
= Plim / / (j^) 2 iA q) (dx,ds)= [ [ {j^fdxds. 

k ^°°Jo J\i^\<^r hi +S JO J\jM\<V6 

Since a.s. Jj 7 ( q )| <1 (7^) 2 <ia; < oo, the last term above tends a.s. to zero as 5 — > 0, 
which implies that Plim Sn) = 0. 

Since then from any subsequence of S n we can extract a sub-subsequence tending to 
zero in probability, we in fact have that the whole sequence Sn^ — > in probability, 
as we need. □ 



Proof of Theorem \4-2\ We decompose ICt,u — ICt into the sum of five terms 
and we show that each term tends a.s. to zero, as n — > oo. We need some further 
notation. Recall that for each q = 1, 2 



and denote 



D[ 9) = [ atf*ds+ [ ai q) dW^\ 
Jo Jo 



so that we have x[ q) = Y t {q) + J$, q = 1, 2. 

Adding and subtracting ^^i^/^'^A/wi^grj^^^'^A^PJj^^r,,} from 
ICt,h — ICt, w e reach 

\ICt,u — ICt\ 
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|^(A i y( 1 )+A i j«)i {(Aj . x(1))2 < r . fe} (A i y( 2 )+A i Jf)i {(A: 



<|5>y«i t(A ,yd))2<9r fe } A ^ (2)l 
i=i 



+ 



Ej=l A i y(1)A i^ (2) ( 1 {(A J X(i)) 2 < r ^l{(A J X( 2 )) 2 <r h }- 1 {(A j y(i)) 2 <9r h } 1 {(A 3 y( 2 )) 2 <9r h ^ 



+ 



{(A J X(D) 2 <r, l } 1 {(A J X( 2 )) 2 <r h } 



+ 



{(A J -X(D) 2 <r4 i {(A j X( 2 )) 2 <r h } 



+ 



(12) 



+ | ^ Aj J^-'Aj J 2 (2) l| (AiX (l) ) 2 <T . h |l| (AiX (2) ) 2< rh } 



The first term tends to zero in probability by Corollary I3.5I The second term 
coincides with 



{(A J -X(i)) 2 <r h ,(A j X( 2 )) 2 <r h ,(A j y(D) 2 >9r h } 



+ 1 



{(A J -X(D) 2 <r h ,(A,X( 2 )) 2 <r h ,(A i y( 2 )) 2 >9r4 



-1 



{(A J X(i)) 2 <r h ,(A J X( 2 )) 2 <r h ,(A j y(i)) 2 >9r h ,(A j y( 2 )) 2 >9r h } 



"l{(A J X( 1 )) 2 >r h ,(A j y( 1 )) 2 <9r h ,(A j y( 2 )) 2 <9r h } — l{(A j X( 2 )) 2 >r h ,(A j y( 1 )) 2 <9r h ,(A j y( 2 )) 2 <9r h } 



+ 1 {(A J X(i)) 2 >rh,(A,X( 2 )) 2 >r h ,(A j y(i)) 2 <9r h ,(A j y( 2 )) 2 <9r h } 



(13) 



All these terms tend a.s. to zero. In fact for the first three ones notice that on 
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{(AjX^) 2 < r^iAjY^) 2 > 9r h } we have JT h > \AjX^\ > \AjY^\ - [A^J^j 
and thus \AjJ^ q) \ > \AjY^\ - ^ > 3^- y/r^ = 2^, so that {(A^-X^) 2 
< r ft , (AjF^) 2 > 9r h } C {(A^X^) 2 < r ft , (A^) 2 > 4r h }, g = 1,2, and thus the 
probability that each one of the first three terms of ( {TBI) is non zero is bounded by 

n 

P { Yl 1 {(A J X( 9 ))2<r h ,(A J J^ ) )2>4r4 ^ } ( 14 ) 

for a suitable g. Now on {(A^JfW) 2 < r h , (AjJ^) 2 > Ar h } we in fact have that 
AjN® ^ 0. Actually, since 



2^T h - \A^\ < |(A,J^) 2 | - |A,F^| < ia^i < 



then 



W2Mo g I + |A^'|>|A^| + |A^'|>|A^)|> 



so 



iAi/ri>v^ l-iir, 



since a.s. for sufficiently small h the quantity 1 — K q y2hlogj^/rh is positive, then 
in fact IAjJ^I > 0, so that AjN^ ^ 0. 



So dH]) is dominated by P{AjN^ ± 0, (A J -J 2 ( ' 2) ) 2 > 4r h } which tends to zero 



as /i — ■> by Lemma |7.3[ part 2). 

As for the last three terms of ( fT3l) note that on {(A^Y^) 2 < 9r^} we have a.s., for 
h small such that AjN^ G {0, 1}, 



AjN^ < \Ajji q) \ = lAjY® - AjD^\ < \AjD®\ + \AjY®\ 



< | A 3 D^ | +3yfr<K q J 2/ilogi + 3V^^0, g = l,2, 
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hence, for small h on {(AjY^) 2 < 9r^} we have AjN^ = 0, j = 1, ..,n. Therefore 
{(AjX^) 2 > r^A/t')) 2 < 9r h } C {(AjD® + Aj J^) 2 > C {(A^^) 2 > 
^■} U {(AjJ^) 2 > g = 1,2; however, by Lemma [7.21 part 1), a.s., for small h, 
1{(AjD(9)) 2 >^4 = 0' thus the last three terms of f|T3|) are dominated by 

for a suitable g. However this last term tends to zero in probability, since 

n 1 n 

i=i i=i 

and Epilog I E?=il {(Aj ^> )2> ^ } ] = 2Mog±£; =1 P{(A,J^) 2 > ^} = 0{^nh). 

We now show that the third and fourth terms of the right hand side of (TT2l) . 
which are similar, tend to zero in probability We have 

YTj=l A i^ (1)A J^2 2)l {(A J X(D)2<r4 1 {(A j X(2))2<r h } 



1 {|A J X(i)|<^JA,jf ) |<2 v ^} 1 {|A j X(2)|<^,|A J jf ) |<2 v ^T} 



+ 1 {|A J X(i)|<^,|A J j( 1) |>2^} 1 {|A J X(2)|< v ^,|A J jf ) |>2 v ^T} ' '° ' 



' -L J' I A . Yd) \<r /^IA.f( 1 )|^0 /ETT.1 



{|A j X(i)|< v ^JA J ^ 1, |>2^} i {|A J X(2)|< v ^,|A,J^|<2 v ^:} 
+ l {\A ] XW\<^,\A j J^ ) \<2^} l {\A ] X^)\<^\A ] jf ) \>2^}_ 

As before on {(AjX^) 2 < r h , (AjJ^) 2 > 4r h } we have that AjN® ^ 0, so, for 
each one of the last three terms of (fl5|) . the probability it is different from zero is 
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dominated by 



Y^PiAjN® ? 0, (A,.J^) 2 > 4r h } -> 0. 



j'=i 



Now we show that the first term of ( Tl5|) is asymptotically negligible. Notice that on 



{\AjX^\ < y/rji, \AjJ\ q) \ < 2^¥l} a.s. for small h we have AN) q) = 0; in fact a.s. 
for small h we have AjN^ e {0, 1}, and 



A 3 -iV (9) < {Ajjf'l = \ AjX {q) - AjD {q) - AjJ. 
< v ^ + sup|A iJ D^|+2V^^0, 



(9)1 



for all j = 1, .., n, for q = 1,2. So we have 



Plim 

n— >oo 



Va-f (1) a-/ (2) i 

Z^^J* lA J J 2 i {|A 7 X< 1 ) 



{\AjXW\<JfH,\A J 4 1) \<2JfH} L {\A J X^)\<^,\A :i 4 2) \<2^} 



< PlimV|A 1 DWA, J 2 (2) |l ) 



{|A i jf ) |<2Vrv:}- 



3=1 



By the Cauchy-Schwarz inequality, last term is dominated by 



Plim 

n^oo A 



E( A ^ (1) ) ; 



\ 



3=1 



{|A J J^|<2 V ^T} 



< »/ / (aW) 2 ^ PlimA/^r =0, 

/q n— >oo 

by Lemma 17.41 

It remains to consider the last term of (fl2l) . which is rewritten as in (fl5|) with 
AjJ^AjJ^ in place of AjY^AjJ^', so that last three terms converge to zero in 
probability as before. As for the first term 



7(2) 



Plim y^AjJ^AjJ^l 



{(A J X(i))2<r h ,(A J J< 1) )2<4r4 1 {(A J X(2))2<r, l ,(A J jf ) )2<4r h }' 



(16) 



3=1 
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we remark that it is bounded in absolute value by 



3=1 



<PlimV^ 1) V 5 '- 2) = ' 

n— >oo 

by Lemma 17.41 p 
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